Equation of motion approach to the solution of Anderson model 
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Based on an equation of motion approach the single impu- 
rity Anderson model(SIAM) is reexamined. Using the clus- 
ter expansions the equations of motion of Green functions 
are transformed into the corresponding equations of motion 
of connected Green functions, which provides a natural and 
uniform truncation scheme. A factor of two missing in the 
Lacroix's approximation for the Kondo temperature is gained 
in the next higher order truncation beyond Lacroix's. A quan- 
titative improvement in the density of states at the Fermi level 
is also obtained. 



The single impurity Anderson model(SIAM)El is one 
of the most fundamental and probably best understood 
models in the field of strongly correlated electronic sys- 
tems. It was proposed to describe the properties of mag- 
netic impurities in non-magnetic metallic hosts (for a re- 
view see e.g. Ref. g). Yet, the SIAM has been used widely 
to imitate mixed valence and heavy fermion systems.ETEl 
Since the SIAM was proposed, a variety of standard tech- 
niques have been applied to it and new methods have 
been developed to study its static and dynamical prop- 
erties, basically in the whole parameter space. ( For a 
recent review see, e.g. Ref. 0.) 

The method of equations of motion(EOM) of Green 
functions(GFs) is one of the most important tools to 
solve the model Haniiltonian problems in condensed mat- 
ter physics. One of its most appealing features is that it 
can work in the whole parameter space. Many authors 
haijEjapplied the EQMj^proach to the SIAM for finite 
[/EIJ3'B or infinite C/.E£rtfl In order to close the chain of 
equations it is usual to introduce the conventional Tyab- 
likov decoupling scheme.Ej Lacroix enajplpyed thus decou- 
pling scheme for higher order GFs.E3'tj In the limit of 
strong intra-atomic Coulomb interaction he did not get 
the correct expression of Kondo temperature i*9in his 
high- and/or intermediate-temperature solutionalB at his 
approximation level. 

In the present paper we shall go beyond Lacroix's ap- 
proximation in a slightly different way. Iij-4he following 
we use the correlation dynamics approachllZl to treat the 
problem. After writting down the hierarchy of EOM of 
the GFs the Tyablikov decoupling scheme is not applied 
directly. Instead systematic cluster expansions are em- 
ployed, which express the higher order GFs in terms of 
same order connected GFs and lower order GFs, the hi- 
erarchy of EOM of the usual GFs is thus transformed 



into that of EOM of the connected GFs. The connected 
GFs are defined such that they can not be reduced to 
the low order ones by any way of decoupling. The hi- 
erarchy of equations of the connected GFs provides a 
natural and uniform truncation scheme. Our formalism, 
which is essentially equivalent to the Tyablikov decou- 
pling scheme, is more sp^stematic. The well known^fesults 
like mean field Lheoryld Hubbard-I appmyiniationjlj two- 
peak solutions ,E3 and Lacroix's resultal3llj are recovered 
exactly in successively higher levels of truncation. To go 
beyond Lacroix we introduce even higher order trunca- 
tion. The Kondo temperature obtained is in agreement 
with the exact onellj except for a prefactor. The density 
of states(DOS) of the SIAM is calculated numerically at 
finite temperature. A quatitative improvement on the 
DOS at the Fermi level is also obtained. p. 

Consider the Hamiltonian of the conventional SIAMtl 

^ = ^ ekfika + ^ £dnda 

ka o" 

+ -^U^hdanda + V^{c'l^da +d'^Cka)- (1) 

o" ka 

The notation used is as usual. In this paper the DOS for 
conduction electrons is taken to be a constant, p{e) = jj^ 
as —D < e < D. The width of the virtual bound state is 
A = 7r{V')/{2D)^ 

Using the Zubarev notationlij Ga.b{^) —^ A; B ^, 
it is straightforward to write down the EOM for the one- 
particle GFs 



V- 



-ed-V] ) <C dy, d+ >== 1 

^ — ' UJ — El. I 



' LO - £k 



and for the high order GFs 



VU -C fldada^d'^ >, 



(2) 



(CJ + Ed - £k' - £k) < d'^Ck'aCka;d^ > = 

-U «: ndadtck'sCkcy\d'l > +y(< d^Ck'g-da;d'^ > 

+ < fidsCkc,;d'^ > - ^ < cl„gCk'sCka]dc >), (3a) 

k" 

{uj + Ek' - £d- £k) < c^'ffrfffCfccr; rfj: >= 

U < nda4,^d^Cka]d+ > +V{-^ cl,gdsdy,dt > 

- < fidsCka\d'^ ^ + X! ^ c+ -Cfe"ffCfc<T; d+ >), (3b) 

k" 
(w + Ek' - Sk- Ed) < cl,gCksda\d'^ >= {cl,^Cka) 



+U < ndsc+-Cksd„;d+ > +V{<€. 4,-d^d„;d+ > 

- < d^Ckgda;d+ > +^ < c^/<jCfcffCfe"^;d+ »). (3c) 
fc" 

To save space, the EOM for the high order GFs ^ 
ndsd„\d+ >, < hdsCka]d+ >, <C dtck9d„\d+ > 
, and ^ c^-dgda'jd^ ^ are omitted. One can find 
their EOM in Ref. 03. Instead of employing directly the 
Tyablikov decoupling scheme, we make use of a cluster 
expansions to separate the connected part of the GFs. 
As an example, the high order GF <^ hdadcr',d^ :s> is 
expressed as follows: 

< ndada,d^ >= Uda <-d^,d^ > 

+ '^hdada,d+ :^c, (4) 

where ^ ... »c represents a connected GF and Uda — 
{hda)- It is straightforward to derive the EOM for the 
connected GFs. We write down the first two as follows: 



£d - Un.ds - ^ 



T/2 



U) - Ek 



■€:dy,d+ >= 1 



+ [/ <C hdgda;d^ >c, (5a) 

[UJ -Ed- C/(l - Uda)] < ndgda;d+ >c 

== Uudai^ - rids) < d„\d'^ > +V^^(< fldsCka^d'^ >c 

+ < 4cfc<Trf<T; d+ >c - <C 4^dada; dt >c). (5b) 

It is not difficult to obtain the EOM of the other con- 
nected GFs such as <C 

ndcrCka;d+ >c,< clgdgdcr;d+ >c, < dtck9da;d+ >c 

, < d'^Ck'ffCkcr; d+ >c, <C cl,gdsCka]d+ >c, and 

In the following we discuss briefly the solutions at dif- 
ferent orders of approximation. The lowest order trun- 
cation, i.e. <^ hdgda'jd^ S>c~ 0, leads to the mean field 
theory, while next higher order truncation, namely, 

< c^gdffdcr;d+ >cW 0, < d^Ckada] d+ >cW 0, and 

<^ ndaCka'id'^ ^c~ 0, corrcsponds to Hubbard-I ap- 
proximation. Subsequently, the solution with two peaks 
localized at £d and £d + U with weights 1 — Uda and 
Uds respectively, can be obtained from the truncation 

< ct^d^dy,d+ »c« 0,< dtckada;d+ >c« 0,< 
dtck'aCka;d+ >c« 0, and < c^dgCkcr;d+ >c« 0. 
To reach Lacroix's approximationjEHTthe connected GFs 
like < dtck'aCka]d+ >c,< c'^,^daCka;d+ >c, and 
<^ c^,g.Ckgdcr',d'^ ^c invovling two conduction electrons 
are neglected. After some algebraic manipulation the GF 
^ drj,di: ;» in the limit of L/ — » oo reads 



<d<,;d+ >= il-Uda-VY^ 



{dtcks) 

LU - £k 



[uj - Ed + «A 



■^ '.1 — ^.- ^^ (,I — Pi., ^ I,! — P7. ■ \ ' 



Eq. (g) recovers exactly Lacroix's results. We shall discuss 
in detail the solution of eq. (^) together with the solutions 
of the next order approximation. 

To go beyond the Lacroix's approximation it is nec- 
essary to consider the EOM of the connected GFs 
< d'^Ck'aCka;d+ >c,< cl,-dsCka\d+ >c, and <C 
c%gCksda]d'^ ;2>c, and to assume higher order correla- 
tion Green's functions to be zero. After a lengthy but 
straightforward calculation, ^ da-;d'^ ^ is obtained fi- 
nally as 



<d^:d+ >= 



1 - nda -VJ2i 



(d-tcks) 



uj-Ed + iA- J2k,k' 



VH<,-Ck-,) 



+ 7 



with 



V"^ s--^ {dtcks) Y^ jdtck'a) 

nda ^ ^ ~ Ek ^ UJ- Ek' 
k k' 



(7) 



(8a) 



-iAS + (3 + 2 



E' 



V^d+Ck^) 



tt' (^ - £k){uj - Ek') 



+ 



T/2 

— E 



{{ndaCt,gCka)c - (dtcka)(4,^d^)) 



UJ ~ Ek 



Xl. "V ^^''^'>^^k"a'^k'a) + (fldadtck's) 
Uda ^T^ ,, 



{uJ - Ek){uJ - Ek') 



ndc 



■E 

k,k' 



{d+dtck'gCka) ~ {dtcks){dtck'a) 



UJ - Ek 



(8b) 



where 



k,k 



^ UJ - Ek ^ UJ - Ek' ■^ UJ - Ek 



g^Y^Y^ {dtdtck'sCka)c - {d+ C+-dgCka) c ,^. 

rida ff, UJ - Ek 

k.k' 

In the above derivation, we have used the Hermitian con- 
ditions like (d'^Cks) — {Ckgdg), etc. To simplify the 
solution and to obtain some analytical results, we no- 
tice the following facts: (i) (d^d'^Ck'g-Cka) ^ and 
{fidcrdg Cka) '^ in the limit of [/ — * cxd, in which dou- 
bly occupied states are unfavorable, (ii) {fidaCk'^Cks) is 
a new correlation function related to {dgfidsCks) by the 
EOM and the spectral theorem. From observation (i) it is 
evident {d^fidsCks) ^ 0. Thus, the spectral theorem im- 
plies that {c^, gfidaCks) ^ 0. The cluster expansion yields 
simply {ndaCk'iyCkff)c-{dtcka){Ck>ffda) '^ -ridsicl'a^^ks)- 
(iii) /3 contains two higher order equal time correlations. 
One is the double hopping correlation {d'^dgCk'gCka)c 
which is approximately equal to — {d'^ Cka) {dt Ck' s) based 
on observation (i). The other is the spin- flip corre- 
lation {d'^c'^,^dgCka)c, which should be calculated self- 
consistently. For simplicity, we just assume this correla- 
tion to be zero. Thus, we have approximately 

^. V y {dtck,){d+ck'a) 

rida f-^, UJ - Ek 



iv) By using the spectral theorem again, it is easy to 
show that J2k' V{dtck'a) = -iAuda- 

Based on the above observations we obtain finaUy 



<d^;d+ >= il-rida-Y^ 



V{4cka) 
UJ - Ek 



:{u;-£, + ^A(1 - S) _ 2^ ■^'^^^'^'^'^''^ 



k,k' 



UJ-Ek 



- E 



y3(4cfcff)(c+,^Cfc'^) 1 



(9) 



'^'^" kk^k" (^ ~ ^k){(^ - £k') 

where 6 is given by eq. (peg). 

The average functions (c'^ig-Cka), {<^k"s'^k'a) and 
{d^Cka) have to be calculated self-consistently. To dis- 
tinguish our results for ^ d„]d^ ^ from Lacroix's we 
use G^{uj + iri) for the Lacroix's solution and G^ (M+irj) 
for our solution. Following Lacroix, eq. (ph readsEHl 



G^iu; + trj) = 

where 

A{lu + irf) - 



1 — Ud — A{uj + if]) 



w - £d + iA + B{uJ + irf) - 2iI^A{w + i-q) ' 

(10) 



TT j UJ ~ LO — ilf 



■n J Lu' — CO ^ irj 

Similarly, the eq.(0) becomes 



(11a) 
(lib) 



G^iuj + i7]) = {l-nd-A-—)[u;-ed + iA{l - — ) 

nd rid 



-2{B - iAA) + —{B- i/\A)] 
nd 



(12) 



For the sake of simplicity, we have considered the non- 
magnetic case, i.e. n^t = f^di = nd which is half of the 
total d-electron number nj. By the spectral theorem: 



rit^ J fiuj')p{T,uj')duj', (13) 

where p(T,uj) = -|lmGd(aj + irj) is the DOS with fi- 
nite temperature and f{uj) ^ cxpKt^-ifj/T+i ] ^^ ^ermi 
distribution function. 

Eqs.(p^ and (|l0| ) with ( p^ constitute two sets of self- 
consistent equations respectively. They can be solved 
numerically. Before performing the numerical calculation 
we simply consider its analytic solutions in the high and 
low temperature limits. 

At high temperature A{ijj+iTii) is a small correction and 
thus an expansion with respect to A is available. After 
some algebra, one obtains the Kondo temperature 



Tk = D exp 



7r(erf - Ep) 



(14) 



for the Lacroix's solution, and 

7r(ed - Ep) 



Tk ^ D exp ■ 



2A 



(15) 



for our solution. Comparing with the exact expression 
obtained by HaldaneEJ 



Tk^iDAY^'cxp 



1/2 ..„ '^i^d - Ep) 



2A 



(16) 



our result is evidently an improvement. It is also helpful 
to check the low temperature result. At low temperature 
limit, as Lacroix noticed, A(lo -\- irf) must be large, espe- 
cially as Lo is near the Fermi level. We consider the zero 
temperature case. In this case, A^ B can be written as 



Ar ■ X ^^*/ N A"" 1 \ijJ — Ef\ 

A{iu + Z7?) = --G*d{uj) {— + In ^ ' 
B(. + ^,) = - -+ln^^ 



(17a) 
(17b) 



Inserting A and B into the Eq.(12), one obtain readily 



G^d{E, 



riN* 



(Ef) 



2iAyTrG^*{EF) + A/7T' 



(18) 



The expression is exactly the same as that obtained by 
Lacroix. Ej 

the following we solve numerically the DOS by 



In 



cally 

the self-consistent equations (^_2|) and I fiJj } in the Kondo 
regime. The following parameters are considered. The 
total number of the d-electron is taken to be 0.9, which 
determines self-consistently the chemical potential. The 
halfwidth D is assumed to be 1, which defines the energy 
scale. The d-electron level Sd is taken to be and the 
width of the virtual bound state A = O.OID. In Fig. |l| 
the DOS for our solution (solid line) is shown at the very 
low temperature T — lO^^A. For comparison the DOS 
for the solution (10) with the above parameters (dash 
line), which is not available in Ref. |l^, is also presented. 
In two cases the lorentzian resonance peaks are slightly 
shifted. In the Kondo regime better behaviors includ- 
ing larger \ed — Ef\ and a more pronounced Kondo peak 
are observed in our solution. A good, test for the DOS 
is provided by the Friedel sum rulejlj which relates the 
DOS at the Fermi level to the occupation Bumber rit at 



T = OK, i.e. p{T^O,Ef) 



2/7rnt 



i^)E 



This value 

is about 31 using the above parameter values. Our result 
of about 25 is in rough agreement with this exact value, 
while agreement is not as good for Lacroix's approxima- 
tion. 

The temperature dependence of the DOS at the Fermi 
level is plotted in Fig. in the regime of the Kondo 
temperature, which is given by eq. (|l5| ) in our calcu- 
lation. It is evident that our result is comparable to 
the nuHierical renormalization group(NRG) result (dotted 
line).E3 The result (dash line) from the Lacroix's approx- 
imation[eq.(10)] with the same parameters above is far 



below the NRG results. Therefore, the higher order cor- 
relation effects can not be neglected in describing the 
system properities quantatively. 

In conclusion, by using the EOM for the connected GFs 
we obtained the solution of the SIAM beyond Lacroix's 
approximation. The missing factor of two in the expres- 
sion of Kondo temperature in Lacroix's approximation 
is regained, which is in agreement with the exact one 
except for a prefactor. The DOS and its temperature 
dependence at the Fermi level are obtained numerically. 
A quantitative improvement on Lacroix's approximation 
is found, which shows that the higher order correlation 
effects are important for describing the system properties 
quantatively. 
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this reference. 
^^ The temperature dependence of the DOS at the Fermi level 
was presented in Ref. hS. But the author used the param- 
eter "^ ^ =0.1 {Eo is the d-electron level), which does 
not describe the Kondo regime; according to Ref. ECJ, his 
parameters describe the mixed valence or empty orbital 
regime. 



FIG. 1. DOS for our solution(solid line) and for the 
Lacroix's approximation (dash line) at T = 10~"'A. The val- 
ues of the parameters are Ud ~ 0.9, Ed ~ 0.0, A — O.OID, 
and D — 1. The d-electron level and the position for the 
Fermi level are indicated by the arrows. 



FIG. 2. Temperature dependence of p{T, Ep). The param- 
eters except the temperature are the same as those given in 
Fig.§ 
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